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The International Space Station (ISS) living areas receive the preponderance of ionizing 
radiation exposure from Galactic Cosmic Rays (GCR) and geomagnetically trapped protons. 
Practically all trapped proton exposure occurs when the ISS passes through the South 
Atlantic Anomaly (SAA) region. The fact that this region is in proximity to a trapping 
“mirror point” indicates that the proton flux is highly directional. The inherent shielding 
provided by the ISS structure is represented by a recently-developed CAD model of the 
current 11-A configuration. Using modeled environment and configuration, trapped proton 
exposures have been analytically estimated at selected target points within the Service and 
Lab Modules. The results indicate that the directional flux may lead to substantially 
different exposure characteristics than the more common analyses that assume an isotropic 
environment. Additionally, predictive capability of the computational procedure should 
allow sensitive validation with corresponding on-board directional dosimeters. 
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Nomenclature 

= magnetic field intensity vector 
= proton kinetic energy 
= distribution function normalization factor 
= altitude 

= upper atmosphere scale height 
= magnetic field dip angle 
= directional proton flux 
= omni-directional (integrated) proton flux 
= parameter defined by Equation (4) 

= Earth radius 
= proton gyroradius 
= parameter defined by Equation (6) 

= pitch angle with respect to magnetic field 
= azimuth angle with respect to magnetic field 
= standard deviation of pitch angle 
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I. Introduction 

The ISS at the present time has evolved as a near-Earth space habitat suitable for continuous human occupation. 
Further evolution of ISS should render it as facility forming a vital part of an expanding space exploration 
infrastructure. This study will look at the radiation exposure aspect of astronaut health and safety by utilizing 
analytical procedures for determining ionizing radiation dose with a view toward implementation as a means of 
shield augmentation for the habitation modules. A CAD model of the ISS 1 1-A configuration specifically dedicated 
to exposure analysis has been developed for this study. 

The first step in the analytical process begins with establishment of an appropriate environment model. For the 
Low Earth Orbit (LEO) environment, the most important contributors to deposition of ionizing radiation energy are 
the trapped protons and the GCR. The present study addresses only the highly directional (vectorial) proton flux, 
which very roughly constitutes about half the total cumulative exposure for long duration missions. However, 
instantaneous dose rates are very much higher during the approximately 10-15 minute SAA transits for which most 
of the trapped proton exposure occurs during a 24-hour day. During the transits, both omni-directional and vector 
proton flux vary from near zero to maximum values, and directionality is controlled by the vehicle orientation with 
respect to the magnetic field vector components. Consequently, an added degree of complexity is introduced with 
the time variation of proton flux spectra along the orbit, for which individual transport properties through the shield 
medium must be taken into account. The deterministic high energy heavy ion transport code HZETRN 1 , developed 
at NASA-Langley, is used to describe the attenuation and interaction of the LEO environment particles along with 
the dosimetric quantities of interest. The ISS geometry defined by the CAD model is finally used to calculate 
exposures at selected target points within the modules, some of which represent locations of thermo-luminescent 
detectors (TLDs). 


II. LEO Environment and Proton Transport 

This section describes the radiation environment selected for the present study and its spatial variation in the 
SAA region. Nominal ISS orbital conditions are prescribed as 400 km altitude at 51.6-degree inclination. Simple 
circular orbit equations have been used to tailor the SAA transits for passage through peak flux regions. Time 
variation of the exposure is defined by these transits. 

A. SAA Protons for ISS Transit 

The standard NASA trapped proton model AP8MIN 2 has been chosen to define a near-worst-case scenario for 
the fluxes. Fig. 1 depicts the orbital tracks in ascent and descent passing through the high flux regions. 



Figure 1. Ascent and descent orbital tracks for ISS through the South Atlantic Anomaly. Symbol spacing 
represents 1-min. intervals; flux contours are in units of protons(>100 MeV)/cm 2 -sec. 
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The differential flux spectra obtained from the environment model are plotted in Figs. 2a and 2b for selected 
points near the region of peak flux. The chosen points are identified by time values in minutes elapsed after 
ascending node point. 




Figure 2. Omni-directional differential flux spectra obtained from the AP8MIN model in central region of 
SAA for (a) descending track and (b) ascending track. 
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The complex low-energy behavior in the proton spectra is not readily explained and is most likely due to several 
influences. Since only higher energy protons (> ~50 MeV) penetrate the ISS structure, the low energy fluctuations 
are unimportant. In order to introduce directionality into the flux spectra, the local magnetic field properties become 
a major factor in the environment. 

Near a mirror point, the spiraling particle paths are nearly normal to the field lines (i. e., pitch angle approaches 
90°). A good account of the theoretical basis for the vector flux of protons in the SAA may be found in Heckman 
and N a kail o’, and computational models have been developed for analyzing the effects of directionality 4,5 . Using 
critical assumptions and approximations, an expression for the directional flux has been found 3 in terms of local 
magnetic field vector, B; altitude, H; ionospheric scale height, h s , and the pitch and azimuth angles (9 and X). This 
formula, in the nomenclature of Kern 4 5 , is expressed as a ratio of the vector flux to the omnidirectional (integrated) 
value: 
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where I is the magnetic dip angle, and r g is the proton gyroradius given by 


sin^v^ 2 

+ 1876£ 

30 

B 



with the proton kinetic energy, E, in MeV and magnetic field strength, B, in gauss. The standard deviation of pitch 
angle is given by 
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with R® representing the earth radius. F N is a normalization factor, parameterized by Kern 5 as: 


F n = (.075/ 8533 + x)exp(-x) 


(5) 


where 


r cos / 

x = — 
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( 6 ) 


When the omni-directional flux is redistributed according to the distribution function of Equation (1), a pattern 
emerges in which most particles are directed in a very pronounced band of zenith and azimuth angles. 


B. Energetic Proton Transport in Shield Medium 

The spectra of Figure 2 have been used as input to the HZETRN code to compute transport through thickness 
ranges of shield material (Al). Subsequent exposures in simulated tissue (H 2 0) are evaluated as dose equivalents 
using ICRP 6 quality factors for normally incident flux on semi-infinite slab geometry. The NASA-Langley 
HZETRN code is a well-established deterministic procedure allowing rapid and accurate solution to the Boltzmann 
transport equation. Details concerning the interaction and attenuation methodology are described at length 
elsewhere 1,7 . Figures 3a and 3b show the resultant dose vs. depth functions obtained from the transport calculations 
that are used to evaluate ultimate exposures at target points within complex shield configurations defined by the 
desired CAD solid model of the full-scale geometric structure. 
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Dose Equivalent Rate, mrem/min 




Figure 3. Dose vs. Depth functions calculated for Aluminum slab geometry at selected times during SAA 
transit: (a) descending track and (b) ascending track. 
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III. CAD Solid Model of ISS 11-A Configuration 

The primary components of the ISS 11-A configuration are the U. S. Destiny Lab Module, the U. S. Unity 
Connections Module (Node 1), the U. S. Airlock and the three U. S. Pressurized Mating Adaptors (PMAs). Also 
included are the Russian Functional Cargo Block (FGB, or Zarya), the Russian Service Module (SM, or Zvezda, the 
Russian Soyuz Spacecraft, the Russian Progress re-supply vehicle, the Russian Docking Compartment and the truss 
structures. A simplified model of this configuration has been constructed for dedicated shield analysis using the 
commercially available CAD software I-DEAS This model consists of 460 separate components, each having its 
own dimensions, orientation, and density distribution defined in near conformity with the actual hardware. A large 
part of the inherent shielding for the astronauts results from the distributed micrometeoroid shield and the pressure 
vessel itself. The cargo in the primary modules also provides additional shielding. In this analysis it is assumed that 
these components are primarily made up of aluminum. A description of a predecessor (configuration 7 -A) of the 
present model may be found in Hugger et al. s Figures 4, 5, and 6 show an external view of the 1 1-A CAD model as 
it appears on a computer screen and split-view illustrations of the 6 target points chosen for this analysis. 



Figure 4. External perspective view of CAD Modeled ISS 11-A configuration. 



Figure 5. Split view of U. S. Lab Module showing selected target points. 
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Figure 6. Depiction of selected target points in Russian Service Module. 

The distributions of thickness for 970 directions has been evaluated in terms of the scaled thickness in g/cm 2 for 
each of the 6 chosen target points for a spherical coordinate system with origin at the point. The ray directions are 
determined for 22 polar angles and 44 azimuth angles plus 2 separate polar angles at top and bottom. The spherical 
coordinate grid is defined so that each directional ray subtends a constant solid angle. The cumulative distributions 
are given for the 6 points in Figure 7. 



Figure 7. Cumulative thickness distribution for selected target points in ISS 11-A configuration. 
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IV. Results of Calculations 


Table I. Calculated Dose Equivalent Rates (mrem/min) for Selected Target Points for Isotropic and 

Directional SAA Proton Environments 

DESCENDING TRACK 

RACK01 LABI LAB4 

Time step Directional Omni Directional Omni Directional Omni 

53 

54 

55 

56 

57 

58 

59 


0.72 

0.61 

0.45 

0.45 

1.09 

0.97 

1.29 

1.11 

0.67 

0.79 

1.99 

1.80 

1.44 

1.29 

0.70 

0.92 

2.41 

2.21 

1.87 

1.78 

0.92 

1.27 

3.63 

3.36 

1.91 

1.95 

1.00 

1.37 

4.01 

3.77 

1.43 

1.58 

0.84 

1.10 

3.40 

3.27 

0.88 

1.07 

0.58 

0.72 

2.36 

2.41 


NODE1 1 


SM5 


SM6 


Time Step 

53 

54 

55 

56 

57 

58 

59 


Directional Omni Directional Omni Directional Omni 


1.44 

0.88 

0.63 

0.63 

1.26 

1.18 

2.47 

1.59 

1.11 

1.16 

2.41 

2.10 

3.07 

1.98 

1.25 

1.33 

2.89 

2.55 

4.66 

3.03 

1.66 

1.80 

4.16 

3.74 

4.79 

3.35 

1.71 

1.94 

4.66 

4.55 

4.05 

2.93 

1.28 

1.54 

3.91 

3.67 

2.90 

2.15 

0.77 

1.01 

2.83 

2.98 


Time Step 

77 

78 

79 

80 
81 
82 


ASCENDING TRACK 


RACK01 

LABI 


LAB4 


Directional Omni 

Directional 

Omni 

Directional 

Omni 


0.67 

0.72 

0.41 

0.49 

1.17 

1.64 

1.10 

1.21 

0.71 

0.82 

1.90 

2.66 

1.41 

1.57 

0.95 

1.08 

2.39 

3.32 

1.39 

1.56 

0.96 

1.07 

2.35 

3.27 

1.10 

1.24 

0.77 

0.86 

1.86 

2.58 

0.63 

0.72 

0.44 

0.49 

1.10 

1.54 


NODE1 1 


SM5 


SM6 


Time Step 

77 

78 

79 

80 
81 
82 


Directional Omni Directional Omni Directional Omni 


1.88 

1.47 

0.72 

0.68 

1.82 

2.15 

3.01 

2.37 

1.19 

1.16 

2.98 

3.20 

3.73 

2.96 

1.55 

1.52 

3.74 

4.13 

3.66 

2.92 

1.51 

1.51 

3.61 

3.81 

2.89 

2.30 

1.20 

1.21 

2.91 

3.01 

1.73 

1.36 

0.68 

0.69 

1.68 

1.97 


8 


American Institute of Aeronautics and Astronautics 






Each of the entries in the preceding table represents the solid-angle integration of dose equivalent rate resulting 
from protons incident on the target point from all directions. Even though the total doses are of the same magnitude 
for both isotropic and vectorial external environments, the directional properties of the radiation field may be vastly 
different for the two cases. This is illustrated in Fig. 8 for the target point designated RACK01 as spherical 
coordinate angle contour plots of the directional dose. 



0 30 60 90 120 150 180 210 240 270 300 330 360 

Azimuth angle, deg. 



0 30 60 90 120 150 180 210 240 270 300 330 360 

Azimuth angle, deg. 

Figure 8. Contour plots of directional dose equivalent as functions of spherical coordinate angles about 
target point RACK01 for isotropic environment (top) and directional environment (bottom). Units are in 

mrem/(min-sr). 
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V. Analysis of Results 

The contour maps of Fig. 8 portray the differences in directional dose distribution and illustrate quantitavely the 
angular variation of exposure intensity. However, such renditions are difficult to interpret and diagnose analytically. 
Present 3-D computer graphic visualization techniques may be implemented to provide displays that lend 
themselves to much more convenient and rapid interpretation. 




Figure 9. Computer-generated distributions of dose equivalent on spherical surfaces centered on target 
point within ISS CAD model for isotropic environment (top) and directional environment (bottom). 
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The illustrations shown in Fig. 9 represent the application of visualization software exhibiting color-coded 
patterns of directional dose mapped onto a spherical surface. The example chosen is a point near that designated as 
SM6 in a relatively lightly-shielded region of the Service Module. The mapping is for a time step on the ascent path 
and demonstrates a case for which the isotropic and directional doses contrast markedly. Such images clearly show 
the impact of the normalized distribution function that results in a re-direction, or “focusing” of the isotropic flux. 
Consequently, in some cases the integrated dose in the directional case may be substantially less than for that of the 
isotropic environment. In other cases, the reverse may occur as may be seen in the tabular results. Such variations 
arise because of the complex interactions of the charged particle environment with the local magnetic field and the 
changing orientation of the vehicle structure. 


VI. Summary and Conclusion 

The primary purpose of this study is to demonstrate by realistic simulation a procedure for accurately analyzing 
and predicting radiation exposures in the confines of a shielded spacecraft. The procedure described can readily be 
implemented in comprehensive specific mission analyses and shield design efforts. In the present study, we have 
attempted only to portray results pertaining to the exposures encountered by ISS in transit through the higher flux 
regions of SAA. A more detailed analysis along these lines would necessarily address the more realistic 2 or 3 SAA 
transits per day of ISS over an extended time period. Near-term plans are to progress from spatial/temporal 
simulation to real-time analyses as directional dosimeter data becomes available from ISS. Such validations will 
provide a stringent test of the adequacy of the theoretical developments and serve to quantify the predictive 
capabilities as they may apply to future human missions as well as to remote sensing platforms. 
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